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Using the finite simulation-cell homogeneous electron gas (HEG) as a model, we investigate the 
convergence of the correlation energy to the complete basis set (CBS) limit in methods utilis- 
ing plane-wave wavefunction expansions. Simple analytic and numerical results from second-order 
M0ller-Plesset theory (MP2) suggest a 1/M decay of the basis-set incompleteness error where M is 
the number of plane waves used in the calculation, allowing for straightforward extrapolation to the 
CBS limit. As we shall show, the choice of basis set truncation when constructing many-electron 
wavefunctions is far from obvious, and here we propose several alternatives based on the momentum 
transfer vector, which greatly improve the rate of convergence. This is demonstrated for a variety of 
wavefunction methods, from MP2 to coupled-cluster doubles theory (CCD) and the random-phase 
approximation plus second-order screened exchange (RPA-I-SOSEX). Finite basis-set energies are 
presented for these methods and compared with exact benchmarks. A transformation can map the 
orbitals of a general solid state system onto the HEG plane wave basis and thereby allow application 
of these methods to more realistic physical problems. 

PACS numbers: 71.10.-w,71.10.Ca, 71.15.-m,71.15.Ap 



I. INTRODUCTION 

The exact wavefunction for the A^-particlc non- 
relativistic electronic Schrodinger equation can be ex- 
pressed as an expansion of Slater determinants which 
span a complete A'^-electron space in which the problem is 
posed. These Slater determinants, in turn, are comprised 
of the antisymmetrized products of spin orbitals, the set 
of which form a complete one-electron space. In general, 
however, neither the complete A'^-clcctron space, nor the 
complete one-particle space can be represented exactly 
and unavoidably we must make do with M spin orbitals, 
and, at most, the corresponding determinants in the 
A^-electron Fock space that these spin orbitals can con- 
struct. 

Even within this finite set of determinants, it is ex- 
traordinarily difficult to construct exact solutions and in 
practice one has to resort to approximate theories which 
in quantum chemistry form the set of standard models^. 
These range from the single Slater determinant used in 
Hartree-Fock theory to the variationally optimised linear 
combination of the full set of Slater determinants found 
by Full Configuration Interaction (FCI}^i^. The coupled- 
cluster and many-body perturbation series form two dis- 
tinct hierarchies. The ground state energy retrieved by 
FCI is the variationally lowest that can be achieved from 
this one-electron basis, within the wavefunction ansatze 
prescribed, and so is often termed the exact solution in 
this basis. 
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However, the true solution to the Schrodinger equation 
can only be reached using FCI in the limit that the fi- 
nite one-particle basis spans all of space which typically 
entails M — ;> oo. Since this limit can never be reached 
in practice, schemes must be devised to find behaviour 
of expectation values to allow for extrapolation to this 
limit. The complete basis set correlation energy, the dif- 
ference between the HF and FCI energies in the limit of 
M — > cxD, is an important goal in ab initio electronic- 
structure theory. Here, we will concentrate on the con- 
vergence of the correlation energy noting that the conver- 
gence of the Hartree-Fock energy and orbitals is generally 
well-understood and in the case of real systems can be 
obviated with pseudopotentials or carefully chosen atom 
centred basis sets. However, the convergence of the cor- 
relation energy in a plane wave basis set, which has sub- 
stantial contributions from electron-electron cusps, has 
not been widely investigated. 

In studies of molecular systems, the CBS correlation 
energy can be reasonably well approximated by extrap- 
olation. In doing so, a certain functional form of the 
correlation energy is assumed, which can be rationalised 
by a partial wave analysis of the wavefunction around the 
electron-electron cusp. Most wavefunction based calcula- 
tions of atoms and molecules employ correlation consis- 
tent Gaussian type orbital (GTO) basis sets, first devel- 
oped by Dunning and coworkers, that show systematic 
behavior for many atoms and molecules^i^, converging as 
where X refers to the cardinal number of the ba- 
sis set^. Since this cardinal number refers to a principal 
expansion, the number of orbitals [M) increases as , 
and this convergence is equivalent to 1/M. 

The application of quantum chemical wavefunction- 
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based methods to the soUd state is a young and emerg- 
ing fieldji"— Even within this body of work, most of the 
approaches have relied on a basis set expansion in peri- 
odic GTOs, where the wealth of knowledge on the con- 
vergence properties of these basis sets is well established 
from decades of studies in molecular calculations. Far less 
work has been undertaken on the convergence of dcter- 
minantal wavefunction expansions in a plane wave basis, 
despite presenting a number of advantages when working 
in the solid state. By specifying a single cutoff parameter, 
an arbitrarily large set of linearly independent and intrin- 
sically periodic basis functions can be produced, which 
require no optimisation, are free of basis-set superposi- 
tion error and well describe the nature of delocalised elec- 
trons, which are particularly difficult for expansion in a 
more localised basis. 

Since wavefunction-based theories will inevitably be 
much more computationally expensive, it is imperative 
to develop methods in which the convergence with re- 
spect to the one-electron basis is as rapid as possible. 
Although complete basis set results using extrapolation 
procedures have been presented for systems in a plane 
wave basisi^"— , a systematic analysis and rigorous justifi- 
cation for these schemes is still lacking. Furthermore, the 
question arises as to whether more efficient basis set trun- 
cations exist within the complete plane wave set, which 
allow for a more reliable extrapolation to the complete 
basis set limit. This paper aims at a rigorous investi- 
gation of different extrapolation methods for the homo- 
geneous electron gas (HEG), which is taken to be the 
archetypal solid state model system, in order to extend 
the practicality of correlated wavefunction expansions in 
plane waves. 

The limiting behaviour of basis set convergence is due 
to the inability of determintantal expansions to describe 
the features of the electron cusp and this is independent 
of the precise parameterisation of the wavefunction aris- 
ing from the underlying method. This allows scaling rela- 
tionships to hold across the whole hierarchy of standard 
models. As such, this paper examines the behaviour of 
basis set incompleteness error in plane waves on the cor- 
relation energy of a finite TV-electron gas by use of second- 
order M0ller-Plesset Theory (MP2), where analysis can 
be directly performed and numerically verified to gain a 
preliminary understanding of this error. This is possible 
due to the MP2 correlation energy of a finite electron gas 
being well-defined in spite of the divergent behaviour of 
this energy at the thermodynamic limit. 

We show that a more natural interpretation of basis 
sets in momentum space can be found that relates to the 
momentum transfer vector. This discussion gives rise to 
a new type of basis set truncation that we can use to 
better eliminate basis set incompleteness error in MP2 
and other theories. We then move away from the electron 
gas as a model system to show how these findings can be 
transferred back to real, solid-state systems. We hope 
that this provides the first thorough analysis of basis set 
incompleteness in calculations where a plane wave basis 



set is used and will allow for extrapolations to the CBS 
limit to be found, both more reliably and more efficiently. 

II. AN ANALYSIS OF PLANE WAVE BASIS 
SET INCOMPLETENESS ERROR 

In this section we will use the archetypal model solid 
state system, the homogeneous electron gas, to better- 
understand basis set incompleteness in plane waves. We 
will introduce the HEG Hamiltonian and show how MP2 
theory can be applied to produce an analytic expression 
for the correlation energy approaching the complete ba- 
sis set limit, which we verify numerically. Although it 
is well-known that the correlation energy arising from 
MP2 theory diverges in the thermodynamic limit due to 
long-wavelength excitations as the band gap closes, the 
qualitative cusp behaviour as inter-electronic distance 
goes to zero is inherently captured by short-wavelength 
excitations^^. As such, using MP2 as a model theory for 
correlation provides a good starting-point for our discus- 
sion of basis set incompleteness error. 

A. Using the electron gas as a model system 

The TV-electron HEG simulation-cell Hamiltonian can 
be written; 

where a and /3 are electron indices and the two-electron 
operator Vap is the Ewald interaction, 

vyi is the Madelung term, which represents contribu- 
tions to the one-particle energy from interactions between 
a point charge and its own images and a neutralising 
background^i"— , and O is the real-space unit cell vol- 
ume. Hartree atomic units (a.u.) are used throughout 
and energies quoted are total correlation energies for the 
system considered unless otherwise stated. 

The one-electron basis set is taken to be plane waves, 

^,(x)^^,(r,a) = yie'''^--<5.,,., (3) 

where the wavevectors kj are chosen to correspond to 
the reciprocal lattice vectors of a real-space cubic cell of 
length L, 

k= — (n,m,0, (4) 

where n,m and I are integers and il = is the real-space 
unit cell volume of a cubic cell. 
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In this basis the HEG Fock matrix, is diagonal, and 
the Hartree-Fock determinant is the normahsed, anti- 
symmetrized product of N plane waves with the lowest 
kinetic energy. 



Do = A [V'i(xi)?/>j(x2)...V'fc(xAr)] 

with the energy, 

_! W , N N 

{Do\H\Do) 



(5) 



47r 



(6) 

where the removal of the q = term in the two-electron 
operator has removed the two-electron Coulomb term, 
corresponding physically to the cancellation of the classi- 
cal interaction between the electrons and the interaction 
between the electrons and the neutralising background. 
The remaining terms in Eq. (|6]) are the kinetic energy, 
the exchange energy and the Madelung energy. 
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(a) Correlation energies retrieved as a function of Ef, 
in scaled units of (^)^ a.u. 
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B. Convergence of finite basis MP2 calculations 

M0ller-Plesset (MP) theory attempts to find the cor- 
relation energy of a system by treating the full electron- 
electron interaction perturbatively within Rayleigh- 
Schrodinger perturbation theory. Taking the zeroth- 
order Hamiltonian as the sum over Fock operators and 
the Hartree-Fock solutions as the zeroth-order wavefunc- 
tions, the first order energy is the Hartree-Fock energy. 
This makes the second-order term (MP2) the leading con- 
tribution to the correlation energy of the problem. 

The MP2 energy can therefore be expressed. 



E 



MP2 — 



E 



\{Di\H'\Do) 
Eo-Ei 



(7) 



where H' is the fluctuation operator defined as the differ- 
ence between the Hamiltonian and the sum over the Fock 
operators. The zeroth-order wavefunctions Dj are the up 
to A^-fold excitations of the Hartree-Fock determinants 
into a complete, typically infinite, basis. Truncating the 
basis set at some M plane waves, these determinants are 

now the O (^) rearrangements of N electrons in M 

spin orbitals. Since H' contains at most two-electron 
operators, only the the O [A^^A/] doubly excited deter- 
minants of Do make a contribution to this energy. The 
single excitations of the reference are zero due to Bril- 
louin's theorem, or since there are no momentum-allowed 
single excitations in the HEG. Finally, the zeroth-order 
energies Ei are sums over the constituent orbital energies 
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(b) Correlation energies retrieved as a function of M ^, where 
M is the number of spin orbitals used. 

FIG. 1: MP2 correlation energies for the 14 electron gas 

_ 3 

at Vs = 5.0 a.u. retrieved function of Ej^ ^ and 
(where M is the number of spin orbitals) tend 
towards a linear relationship as the complete basis set 
limit is approached. In each plot, the dotted lines refers 
to CBS limits for each basis set size, which are obtained 
by a linear extrapolation of this point and the previous 

three points (sometimes not visible on the graph). 
Using these extrapolated estimates, it can be seen that 
the power-law is smoother due to fewer finite size 

effects, this is due to M being a more appropriate 
variable to consider how much correlation energy the 

basis set retrieves. For this system, M = 1030 
corresponds to a kinetic energy cutoff of 1.3077 a.u. or 
35.59 eV, which changes with both N and for the 
HEG. 



J Gocc 



4 



where in these equations, i refers to any occupied orbital 
and a refers to the virtual orbitals. These equations in- 
clude an exchange energy explicitly, and in the thermo- 
dynamic limit tend towards the well-known forni^i 



where x = k/kp and 



1 



2x 



In 



1 + x 



1-x 



This allows the MP2 energy to be re-written as, 
1 ^ \{ij\vi2\ab) ~ {ij\vi2\ba)\'^ 



Emp2 - 4 XI 

ij ^occ 
afc^virt 



Cj + fi 



(9) 



(10) 



(11) 



where indices i,j,a and b are spin orbitals. 

This can be solved directly for the non-interacting ref- 
erence in the limit of both an infinite number of electrons 
and an infinite virtual k-space^. However, the limit of 
a finite number of electrons is dependent on the form of 
potential Vi2 and the shape of the real space unit cell. 
Furthermore, it is typical to use a finite basis set to de- 
scribe the virtual manifold, which can be achieved in the 
plane-wave basis with a choice of kinetic energy cutoff, 
Ek ~ \k1 such that, 

E-EE ;E-EE 



zjGocc (7i (7j 0<ki<kf 
0<kj<kf 



ab^virt 



CTi o-j kf<ka<kc 



(12) 

where ki = |ki| etc. and sums over spins have been writ- 
ten exphcitly. Using this substitution, Eq. ([TT|) can be 
re-cast in a finite basis. 



Emp2 - Y1 E 

0<ki<kf kf<ka<ka 
Q<kj<kf kf<kt<kc 

- E E 



2|(kikj|{)i2|kQkf, 



e(kO + e(k,)-e(kJ~e(kb) 



{kjkj |i)i2 Ikgkb) (kjkj |{)i2 Ikfckg) 

(13) 

where the sums over spins have been taken leaving a spin- 
free expression. The two sets of terms are referred to as 
direct and exchange-like terms respectively. 

Defining Vq as the q Fourier component of the poten- 
tial, the four index integrals can be evaluated, 

(14) 



(kjkj|{)i2|kakt) = Wk.-k„<5ki-k„,k,-kt, 



yielding, 



EmP2 = Y1 E 



Jki-k^.kj-k, 



2 <-k„ 



0<ki<kf kf<ka<k^ 
0<kj<kf kf<kb<kc 



■EE' 

0<ki<kf kf<ka<ka 
0<kj<kf k}<kh<kc 



3ki-k„,kj-k(,- 



?^ki-k„i^kj-k„ 



(15) 



where Ae is the difference between eigenvalues and de- 
pends on the four indices. Values of k^ in this represen- 
tation are constrained to obey momentum conservation. 



kj — ka -\- kfo. 



(16) 



due to (5ki-k„,kj-kb, and, therefore, the sum over kf, 
makes at most one contribution for every k^, kj and k^. 
This is the formulation of the MP2 energy that we will 
refer to as the iJfc-cutoff scheme. 

The question we now seek to address is: how does the 
correlation energy captured by MP2 increase with the 
energy cutoff of the basis set on approach to the com- 
plete basis set limit? Since this question has not been 
addressed for plane wave basis sets, it is appropriate to 
conduct a simple analysis as follows. 

We seek an expression for the error of a finite calcula- 
tion conducted at a kinetic energy cutoff Ek 



1 7,2 

2'^c 



A£;MP2(fcc) = -EmP2(oo) - £^MP2(fcc), 



(17) 



where the -EMP2(fcc) is the finite-basis MP2 energy given 
in Eq. (fT5|) . This can be evaluated by changing the limits 
on the sums, such that: 



Ai?MP2(fcc) — X X (^ki-k^.kj- 



Q<ki<kf ka>k^ 
0<kj<kf kb>ka 



2 t-g,-k„ 
' Ae 



- E E^. 

0<ki<kf ka>ka 
0<kj<kf kb>k^ 



k^ — ka ,k^ — k), 



^ki -ka«kj-k„ 



(18) 



It is possible to simplify this expression in the high 
basis set limit. The orbital energies become dominated 
by high-energy kinetic energy contributions, whereupon 
Ae oc fc^. As fca 3> ki and ka 3> kj, the numerator tends 
towards a behaviour of 1 /fc^. In this limit, the summation 
of ki and k^ yields a constant factor, and the Kronecker 
delta reduces the double-sum over virtual orbitals to a 
single sum. 

This leads to a leading-order expression of. 



E 

ka>k. 



AEMP2{kc) PC ^ 



(19) 



where the sum can replaced by a spherically symmetric 
integral and evaluated as. 



A£;MP2(fcc) oc / dfca ^ kl 



cx 



1 

kl' 



(20) 



This is equivalent to E^^ ^ , due to the definition that 
Ek = \k^i or M~^, where M is the number of k-points 
contained within the sphere defined by k^. In passing, 
we note that this is the behavior that is also found 
for the correction to the energy in the random phase 
approximation^^. 
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Figure [T] shows numerical verification of tliis relation- 
ship using rg = 5.0 a.u., a typical of real materials. 

_ 3 

In Fig. [Tal a relatively rapid tendency to follow & E^'^ 
power-law is found. Extrapolated results at each basis 
set (using this basis set size and the previous three ba- 
sis set sizes) show rapid convergence to the infinite ba- 
sis set result, although this tendency is not smooth due 
to shell-filling effects. When instead a power-law 
extrapolation is used, as in Fig. Ilbl this convergence is 
somewhat smoother and better behaved for small basis 
set sizes (when the difference between the two power-laws 
is more pronounced). 



III. MOMENTUM TRANSFER VECTOR 
CUTOFF SCHEMES 

In this section, we develop a different type of basis set 
truncation for the HEG, based on the momentum trans- 
fer vector. Rather than the conventional definition of a 
single basis set for the whole calculation, we take the view 
that the basis set can be defined differently for each elec- 
tron or each electron pair. This definition is not unique, 
even given a spherical cutoff, and we develop three types 
of basis set truncations showing that there is one that 
gives more rapid convergence to the CBS limit. This has 
the physical equivalence in reciprocal space of smearing 
out the rigid spherical cutoff into the surrounding space. 
This is motivated by a physical picture that electron co- 
alescences should be treated on the same footing in mo- 
mentum space. 



A. Introducing the momentum transfer vector 

Considering a general same-spin electron-electron- 
hole-hole excitation ij ab connected by a matrix ele- 
ment (kikj |-Ci2 |kak(,) — (kikj 11)12 Ikfeka) . We can therefore 
define two momentum transfer vectors for the excitation, 
g and g'. 



ka = k, + g 



kfc = ko 



ka = kj-g' ; kb^ki+g', 



(21) 



(22) 



where the allowed g vectors are such that k^ and kf, are 
both not in the occupied manifold. It is possible to re- 
write the sum over k^ and k(, in Eq. (jl5p in terms of these 
vectors, 



0<ki<kf /c/<|ki+g|</cc 
0<kj<kf kf<\-k,+s'\<ka 



(23) 

where similar to before g' is specified uniquely by g, k^ 
and kj using <^g,g'+ki-kj • 



By analogy with previous work in solid-state 
systemsiiii^, we now consider cutoffs that limit the ex- 
tent of the momentum transfer vectors, and as such we 
impose a cutoff on the g-vectors such that they do not 
exceed a kinetic energy Eg = ^g^, and such that ka and 
kh never reach the iS^-cutoff value kc- The upper limit 
in the sum becomes entirely determined by gc'. 



(2 vl - VgVg') 



(2 vl - VgVg,) 
Ae 



EmP2 (gc) - ^ ^ ^g,g'+k.-kj 

0<ki<kf kf<\ki+g\ 

0<kj<kf fcj:<|ki+g' 

(24) 

This gives us a new form of basis set truncation whose 
behaviour in the large g limit might be different to 
Emp2 (kc) (Eq. ([231)), which we will now investigate. 

It is also possible to remove the upper limits on the 
sums, replacing them in with radially symmetric step 
functions in k-space. 



e (g - gc) 



1, |g|<5c 
0, otherwise 



(25) 



yielding, 

EmP2 ~ 



2^ Og,g'+k,-k,- 



0<ki<kf fcj><|ki+g| 
0<kj<kf kf<\ki+g'\ 

^g,g'+ki-k. 

0<ki<kf fej<|ki-|-g 
0<kj<kf fc^<|ki-|-g'l 



Ae 



VgVg'Q{g - gc)Q{g' - gc) 



(26) 



This cutoff is shown diagrammatically in Fig. [2a] for 
a specific electron pair, k^ and kj, illustrating that the 
basis set used to represent the virtual manifold is now no 
longer consistent between different electron pairs. By al- 
lowing {g} and {g'} to span a certain range in reciprocal 
space, the virtual space represented by the sets {ki + g} 
and {kj + g'} span a range dependent on k^ and k^ re- 
spectively. 

More severely than this, the basis set we have de- 
fined is also different for each electron. Considering 
a specific single same-spin excitation (kimi')i2|kakf,) — 
(kikj|?)i2|kbka) this means that sometimes the shorter 
momentum transfer vector is allowed while the longer is 
disallowed. This is shown diagrammatically in Fig. I2bl 

Since this cutoff takes the view that each electron has 
its own basis set, this will be termed the local Eg-cntoS. 
In Eq. pSj) this is represented by the different ranges 
of the sums over the direct and exchange-like terms. In 
the exchange-like term the product of the step-functions 
serves to disallow some longer-momentum events. The 
implication of this is that in the general hole pair func- 
tion space |kakfc) can be allowed while its permutation 
Pi2|kak;,) = — |kf,ka) Can be absent. This implies that 



(a) This is a diagram of the local Eg-cntoS. Using a 
simple momentum transfer cutoff scheme (Eq. | |26| |'1 

makes the virtual space accessible to an excitation 
kijkj — >■ {ki + g},{kj + g'} dependent on ki and kj 
when they arc not at the F-point. This implies that 

the sets {k; + g} and {kj + g'} arc not the same. 



(b) This is a diagram of the local Eg-cutoS. 
Comparison between a specific excitation ki,kj — >■ 
ka,kj and ki,kj ki,,ka. These differ only in the 
permutation of the hole states (or, equivalently, the 
electron states). In the case of the local Eg-cutoS 
cutoff scheme (Eq. II26I I). one term (solid line) is 
allowed and the other term (dashed line) is disallowed. 




(c) This is an illustration of the intersection Eg-cutoff. 
One solution to the problem illustrated in Fig. I2bl is to 
only allow excitations to the region of k-spacc formed 
by the overlap of the two regions in Fig. I2al Now both 
terms are either disallowed (as shown here) or allowed 
and pcrmutational symmetry is restored. 




second and different solution is to extend the allowed 
space of Fig. l2al to anywhere that either {ki + g} or 
{kj + g'} would be allowed, also restoring the 
pcrmutational symmetry 



FIG. 2: Discussion and diagrams of cutoffs using monientuni transfer vectors. The wliite circle represents the Fermi 
sphere, and the volume excluded from the virtual space by occupation effects is not considered. 



not all terms accounted for in the direct term arc prop- 
erly balanced by the exchange-like term, and the anti- 
symmetry of the wave function is ultimately not properly 
restored. 

In the conventional basis set scheme described in 
Sec. Ill B1 there is a variational principle: for a finite basis 
set you are guaranteed to not retrieve more correlation 
energy than the complete basis set limit. As the basis set 



is enlarged, the correlation energy is systematically low- 
ered to the complete basis set limit correlation energy. 
This variationality is broken by use of a local Eg-cutoS. 

We can define two further ways of defining a Eg-cntoS, 
which do not suffer from these limitations. In the inter- 
section £^g-cutoff, we force the direct term to be removed 
from the sum if the exchange-like is rejected for the same 
{g,g'} pair: 
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MP2 - ^ ^ f^g.g'+k. 

0<fci<fc/ fcj<|ki+g 
0<kj<kf fc^<|k,+g'| 



2vlP (g,gO 



E E ' 

0<ki<kf fcy<|ki+g| 
0<kj<kf fe^<|k,+g'| 



Sg-g'+ki-kj ■ 



Ae 
Ae 



where P (g, g') is given by, 

p(g,g')-e(g-5c)e(5'~<?c). 



(27) 



(28) 



which can be thought of as a masking function in that it 
disahows certain electron pairs from being connected to 
different parts of the virtual space. 

In the union iSg-cutoff, we force the exchange- like term 
to be preserved if the exchange-like term is rejected for 
the same {g;g'} pair by use of, 

P (g, g') = e (g - gc)+e {g' - g,)-Q [g - g,) 9 [g' - g,) , 

(29) 

where the term O (g — g^) 6 {g' — gc) prevents double- 
counting when Q [g — g^) and O {g' — g^) are both 1. 

The cutoffs are named after how they are generated 
from the sets {k^ + g} and {kj -1- g'}, shown in Fig. [2c] 
and Fig. [2dl 
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FIG. 3: Comparison of correlation energy retrieved as a 
function of basis set size for a variety of cutoff schemes. 



B. Comparison of the different cutoffs 

Figure [3] shows correlation energies with these different 
cutoff schemes. In the case of the E'g-cutoffs, the M is 
the number of spin orbital basis functions in a sphere 
with radius gc centred at the F-point. In some cases, 
in particular the union i?g-cutofF, the number of basis 
functions used for the calculation is higher, but we believe 
that there is no better parameterization of the size of the 
basis set than this effective M . 

All of the schemes regardless of the cutoff scheme con- 
verge to the CBS limit ultimately as 1/M. As the size 
of the basis set goes to infinite extent, all i^^-cutoffs ul- 
timately tend back towards the Ek picture since the dis- 
placement of the occupied fc-points from the F-point be- 
comes negligible and the lines become identical. For the 
intersection and local £^g-cutoff schemes the curves, how- 
ever, merge only at very large basis sets. 

The positioning of the curves of each Eg basis set can 
now be compared with that of the corresponding Ek ba- 
sis set. Both the intersection and union i5g-cutoffs can be 
thought of as lying in a larger Ek-cntoS. basis set and are 
variational upper bounds of this larger basis set energy. 
In the intersection iJ^-cutoffs scheme, terms are effec- 
tively removed from electron pairs that are of significant 
distance from the F-point. Furthermore, all excitations 
lie within gc of the F-point, meaning that this basis set 
produces a variational upper-bound to the iJ^-cutoff ba- 
sis set of the same size {kc = gc)- In contrast, the union 



iJg-cutoff augments the basis set for those electron pairs 
that are not at the F-point by including basis functions 

that can have as high an energy as ^ {^g^ -\- k^^ . As such, 

this is now a variational upper bound of the Ek basis set 

that completely encloses the radius (^g^ + k^^^ . 

In contrast, the local Eg-cntoS is neither variationally 
bounded by the complete basis set limit nor any Ek-cutoS 
basis set. In general, it can be considered that it has 
fewer exchange-like terms than the corresponding union 
basis set, and as such will produce a lower correlation 
energy than all of the basis sets with the same cutoffs. 
Since the exchange-like terms in the correlation energy 
are positive and partly neglected, the correlation energy 
becomes more negative than for the corresponding union 
basis set. Although this seems advantageous in the first 
instance, as it seems to retrieve a greater fraction of the 
CBS correlation energy, already in Fig. [3] it can be seen 
that there is a tendency for this curve to arc at low basis 
sets, and could even have a maximum point in extreme 
cases. 

Each cutoff has a separate behavior when a 1 /M be- 
haviour is used to extrapolate the result from a series 
of finite basis calculations. As noted previously, the Ek- 
cutoff basis set suffers from strong finite size effects, caus- 
ing the extrapolation to behave jaggedly around the CBS 
result. This can be thought of being due to trying to 
recreate a spherical cutoff with a cubic grid. In common 
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(a) (b) 

FIG. 4: Graphs comparing (a) union Eg-cutoS and (b) Ek-cutoS for a range of different densities. As the density is 
lowered, the extrapolations become more distant from the CBS limit even for MP2 theory, due to a rising 
contribution from exchange in the Hartree-Fock orbital energies that is not well-behaved with respect to M . 
Furthermore, more pronounced finite size effects are seen. The extrapolated results shown by the dotted lines are 
only represented with error bars in two cases (r^ = 0.5 a.u., = 20.0 a.u.) for clarity. Complete basis set limit 
energies from which these basis set incompleteness errors are derived are tabulated in Tabic ID 



with this, the intersection iS^-cutoff has even stronger 
shell-filling effects, which are more pronounced because 
this basis set is trying to recreate the overlap between 
two spheres with this cubic grid, a shape with an even 
smaller volume to surface ratio. The local and union Eg- 
cutoffs have much smoother convergences with 1/M and 
their extrapolated results converge much more smoothly 
to the CBS limit. This could be because we have replaced 
spheres in fc-space with more complex objects, and also 
are summing in more excitations for a given M . 

In conclusion, the union i?g-cutoff seems to have the 
most desirable properties: variationality, correct symme- 
try of the wavefunction, smoothness and speed of con- 
vergence and extrapolation. When the density, as rep- 
resented by Ts, is changed we might expect these rela- 
tionships between the cutoffs to change. In Fig. HI we 



have considered the fraction of the CBS correlation en- 
ergy obtained by basis sets at both higher and lower den- 
sities (rg — 0.5 — 20.0 a.u.) for the union Eg-cutoS and 
Sfe-cutoff. As rg is raised, the basis set extrapolation be- 
comes increasingly distant from the CBS result at smaller 
basis set sizes. This is due to the rise in the contribution 
from the exchange-like term in the MP2 energy, which 
has a less well-defined convergence with respect to the 
M parameter that we are using. Furthermore, finite size 
effects become more visible. From these graphs, it is pos- 
sible to see that for this system the union Eg-cutoS does 
continue to be the cutoff of choice for the reasons outlined 
above. For completeness, values of the complete basis set 
correlation energy for MP2 are presented in Table HI 
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Ts (a.u.) 


Correlation energy (a.u.) 


0.5 


-0.575442(1) 


1.0 


-0.499338(2) 


2.0 


-0.398948(2) 


5.0 


-0.255664(4) 


10.0 


-0.163951(6) 


20.0 


-0.09749(1) 



TABLE I: Values of the complete basis set limit MP2 
correlation energy obtained by extrapolation for 
TV = 14. The error estimate, shown in brackets refers to 
the random error in the last digit from extrapolation. 



A. Direct extrapolation of CCD and RPA+SOSEX 

In both CCD and RPA-hSOSEX, the energy can be 
written in a configuration-space formalism as, 



M 



Ecorr (A/) 



(M) 



(30) 



where Xk^' ^^'^ k-point labelled contributions to the 
energy which express a product of amplitudes on double 
excitations of the reference determinant and the corre- 
sponding Hamiltonian matrix element. 



IV. GENERALISATION TO OTHER 
SINGLE-REFERENCE QUANTUM CHEMICAL 
METHODS 



In this section, we seek to generalise the discussion 
above to other single-reference quantum chemical meth- 
ods, in particular the coupled-cluster doubles (CCD) 
theory^! and the random phase approximation plus 
second-order screened exchange (RPA-fSOSEX)^. In 
these methods, the energy estimator depends on ampli- 
tudes, all of which vary when the basis set size is changed. 
This is in marked contrast with MP2 theory, in which 
only those basis functions added when a basis set is en- 
larged acquire new contributions to the energy. We will 
discuss two possible strategies for extrapolating the en- 
ergy to the CBS limit. We will show that a direct extrap- 
olation based on calculations at different basis set sizes is 
one method for achieving the CBS limit, and that both 
CCD and RPA+SOSEX correlation energies behave as 
1/M in common with the MP 2 correlation energy. How- 
ever, this suffers from the same slow convergence of the 
i?fc-cutoff strategy outlined from MP2. 

In an attempt to emulate the more effective extrapo- 
lation to the CBS limit provided by the momentum vec- 
tor transfer cutoff schemes, we introduce a new approach 
to this problem, single-point extrapolation, in which the 
contributions to the energy from a single calculation are 
re-grouped according to their arrangement in recipro- 
cal space to form energy estimates from effective ba- 
sis set sizes. These smaller effective basis set energies 
are then used to provide an extrapolation to the CBS 
limit. Although it will be demonstrated that this ap- 
proach does provide more effective convergence, ampli- 
tude 'relaxation' as the basis set size increases causes a 
problem with CCD at = 5.0 a.u. However, this cru- 
cially also allows for the adaptation of extrapolation to 
solid state systems, where direct extrapolation is not only 
slower to converge but also made difficult to achieve by 
the PAW approximation. 



Xk^;-(2^^k.-k„-«k,-kjt^''. (31) 

Appendix 1X1 contains a more detailled discussion of CCD 
and RPA-f-SOSEX and explains the evaluation of the re- 
spective tfj amplitudes. The additional complexity com- 
pared to MP2 is that all tfj vary when the basis set size 
is increased. Again we quantify the size of the basis set 
of the virtual orbitals using M. M corresponds the total 
number of orbitals inside the cutoff sphere (iJ^-cutoff in 
MP2). 

Performing calculations at different cutoffs allows di- 
rect extrapolation of the energy, with this behaving as 
1/Af in the large M limit. Comparison between the 
finite basis-set energies retrieved by different quantum 
chemical methods is shown in Fig. [S] for the 14-electron 
problem at Vg = 1.0 and rg — 5.0. All methods con- 
sidered show a 1/M relationship in the high M regime. 
In this graph, exact benchmarks from a new electronic 
structure method called initiator full configuration inter- 
action quantum Monte Carlo (j-FCIQMC) are presented 
for comparison from Ref. |2^. This method utilises a 
stochastic algorithm to calculate FCI accuracy energies 
at greatly reduced computational cos t'^^i'^^ . The gradi- 
ent and onset of the 1/M behaviour varies with and 
method. At the higher r^-value CCD best resembles the 
FCI behaviour, with MP2 and RPA-HSOSEX resembling 
one another. All of the methods behave similarly with M 
at the lower Tg-value. The ability of RPA-I-SOSEX to re- 
trieve most of the FCI correlation energy at = 5.0 can 
be attributed in part to capturing too much (dynamic) 
correlation energy at high M. The cross-over between 
CCD and RPA+SOSEX (Fig. [5b| highlights the diffi- 
culties of comparing methods at a finite basis set size. 
Clearly RPA+SOSEX or MP2 can not be used to esti- 
mate the finite size and basis set incompleteness error of 
i-FCIQMC or CCD, whereas CCD may be weU suited 
to correct these errors in i-FCIQMC. No attempt has 
been made to extrapolate these methods to the thermo- 
dynamic limit, in which MP2 is well-known to diverge, 
since this is beyond the scope of this paper. 
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FIG. 5: Comparison of correlation energy for the iV = 14 system for (a) rg = 1.0 (b) Vg = 5.0 retrieved as a function 

of basis set size for a variety of quantum chemical methods. 



B. Single-point extrapolation of CCD and 
RPA+SOSEX 

We now seek a momentum transfer vector cutoff 
scheme for CCD and RPA+SOSEX, in particular aim- 
ing to re-produce the properties of the union iJg-cutoff 
explored in Sec. lIIII After performing a single calculation 
in a basis set, which is now labelled M' for reasons that 
will become clear shortly, 

occ A/' 

i?_,,(M') = ^^xl::^;(A/')- (32) 

ij ah 

Applying the masking function P (g, g') defined previ- 
ously for the union iJg-cutoff, 

P (g, g'; M) - 6 (.g - g,) + 6 (g' - g,) 

-Qi9-9c)e{g'~gc), ^ ' 



which is associated with a basis set size M (described in 
SecHH), to Eq. ^ yields, 

OCC A'l' 

i?co„,off (M', M) = ^ ^ x;::^^^ (mo 

ij ah ^ ' 

X Pg (ki — ka , kj — kfi ; A/) , 

where we have explicitly noted that this formulation of 
the correlation energy is dependent on both M' and M . 
These correlation energies are labelled both by a true 
basis set size M' and what we will call an effective basis 
set size M . We now follow the procedure of performing 
a single calculation with M' spin orbitals, take the am- 
plitudes and apply the relationship given in Eq. (|34p for 
different values of M. 

Analyzing Eq. ([M)) . it is possible to see that there are 
two limiting values for i?corr,eff. When Af = 0, the effec- 
tive basis set correlation energy is zero, and when M is 
such that all possible momentum transfers are included 
in the sum (when g^ > k^ + kf), the effective basis set cor- 
relation energy is simply the basis set correlation energy 

E',,„ (Eq. m)- 
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In between these limits, if the amphtudes tfj 
are always the opposite sign to the matrix element 
(2uki-ka ~ i^kj-ka), there will be a monotonic decrease 
of i?corr,cff {M' , M) to the basis set correlation energy as 
M is increased. 

In MP2 theory, this monotonic decrease will be strictly 
observed, and can be shown to be identical to the union 
-Eg-cutoff scheme when kc > gc + kj. For this region, 
< gc < — kf , therefore the same tendency to follow a 
1/M behavior will be seen. When gc > kc — kf, deviation 
from this behavior will be seen due to momentum transfer 
vectors being disallowed from not being in the original fc^ 
basis. 

Unlike the previous formulation, this can now be ap- 
plied to any method with an estimator of the form 
Eq. ([3T|) . However, since the amplitudes also depend on 
AI' , this is an approximation and convergence with this 
second cutoff should also be obtained. 

Figure [B] shows that these effective basis set energies 
have the property that they also converge as 1/M and can 
be used to extrapolate for a CBS estimate. Fig. [7] shows 
these extrapolations for RPA+SOSEX and CCD, com- 
paring them with conventional direct extrapolation. In 
general, the RPA-I-SOSEX correlation energy converges 
faster using the i?g-cutoff single-point extrapolation than 
the i?/j-cutoff direct extrapolation, which also has the ad- 
vantage that only one calculation needs to be performed 
at a single basis set size. For CCD, this advantage is 
greatly obscured by finite size effects (which would be- 
come less for larger system sizes) and is not seen at all 
for Ts = 5.0 a.u. due to flattening off of the finite basis 
set correlation energies and greater coefficient relaxation 
effects arising from stronger correlation. 

Extensive discussion and analysis of relaxation effects 
are beyond the scope of this paper, but this method has 
also been successfully applied to the stochastic quantum 
chemical method i-FCIQMC, and the further benefits of 
applying such a technique in a stochastic framework are 
discussed in Ref. HI. 




FIG. 6: Comparison between direct extrapolation and 
single-point extrapolation (SPE) for RPA+SOSEX on 
the TV = 14, rs ~ 1.0 a.u. gas. In the conventional 
direct extrapolation, calculations are performed at a 
series of basis set sizes M and then extrapolated using a 
1/M fit to the high M fimit. In the SPE, a single 
calculation is performed at an overall basis set size of 

M', in this case M' = 682, and effective basis set 
energies are constructed according to Eq. ((34)) over the 
full range of M. Some of these points are discarded as 
A/ approaches M' since not all momentum transfer 
vectors can be accommodated within the basis set 
(dashed green line, discussed in the text). The 
extrapolations are shown by dotted lines, and agree in 
the CBS limit within reasonable extrapolation error 
estimates (~ 2 x 10~'^a.u.). 

fore examine the improved extrapolation strategics based 
on the single-point extrapolation scheme discussed previ- 
ously (union Eg-cutoS and i?fc-cutoff), which we believe 
to restore variationality and correct symmetry properties. 
Finally, we show how this can be applied to an example 
periodic system. 



V. APPLICATION TO GENERAL SOLID STATE 
SYSTEMS 

In this final methodological section, we discuss extrap- 
olation schemes available to solid state calculations using 
a plane wave basis set. We start by noting that for solid 
state systems the previous methodology of an £fe-cutoff 
fSec. lIIB)) . or equivalently, an M 'true' basis set fSec lIVj) 
is not easily defined. Previous work to resolve this in a 
plane- wave basis set has used a resolution of the identity 
basis set to identify Hamiltonian matrix elements. Fol- 
lowing a similar argument made in this paper, previous 
authors have found that the correlation energy converges 
with respect to this auxiliary basis set as 1/M. How- 
ever, this greatly resembles the local iJ^-cutoff described 
in this Sec. IIIII and with the most severe penalty being 
that it is not variational with the CBS limit We there- 



A. Formulation of the single point extrapolation 
for solid state systems 

The correlation energy expression in general wavefunc- 
tion based methods is given by 

M' 

£;corr(M') = EE*^(2«^ - ^^)*- (35) 

ij ah 

The indices j, k and a, 5, c, d refer to occupied and un- 
occupied orbitals, respectively and are understood to be 
a shorthand for the band index and Bloch wavevector. In 
contrast to the homogeneous electron gas, however, the 
orbitals are no longer constituted by plane waves and 
correspond to eigenfunctions of the respective Hartree- 
Fock (HF) or Kohn-Sham (KS) one-electron Hamiltoni- 
ans. M' corresponds to the number of basis functions 
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FIG. 7: Correlation energies calculated for the 14-electron system using RPA+SOSEX at (a) — 1.0 and (b) 
Ts = 5.0 and using CCD (c) Ts = 1.0 and (d) Ts = 5.0. Direct extrapolation and single-point extrapolation (SPE) are 
compared at a variety of basis set sizes (SPE curves just show extrapolated results). SPE performs significantly 
better for RPA, due to reduced coefficient relaxation, and at lower correlation strengths. 



used in the description of occupied as well as unoccu- HF/KS orbital. We note that 
pied orbitals. vf^ and V^^ refer to electron repulsion in- 
tegrals and many-electron wavefunction amplitudes, re- 
spectively. 



(38) 



,,a6 



(7^,|r)(r|^,)(^,|r')(r'|V'6) 



dr'dr. 



r — r 



(36) 



For the sake of brevity, we have neglected single exci- 
tation (SE) contributions to the correlation energy in 
Eq. p5p . Depending on the approximation and refer- 
ence determinant used in calculating the wavefunction, 
SE contributions might have to be included but do not 
modify any of the conclusions drawn below. 

We now seek to apply the previously outlined union 
Eg and Ek single-point extrapolation scheme to general 
solid state systems. To this end we introduce a projection 
matrix that transforms the HF/KS-orbitals onto a plane- 
wave basis set and reads 



UnG = (0„|G), (37) 

where |G) is a plane wave, e**^'', and 0„ constitutes a 



If the employed finite plane-wave basis set is complete 
and large enough to span the space of all orbitals 0„, 
UnG becomes a unitary matrix. However, in our case, we 
use fewer orbitals than plane waves. We typically choose 
a plane- wave basis set for which UnG is not full rank, and 
calculate U^^ using a singular value decomposition. 



Inserting Eq. ^E^l into Eq. ([35]) gives 



t: 



,-,G'Gn 



V, 



G'' 



(39) 



G,G',G" 



where 



and 



V^g"^'' =J2U"'G"jU^'aU^},V- 
jab 



(40) 



tiG^' ~ U*ci',UaGUbG'tf!i 
jab 



(41) 
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In contrast to Eq. ((35)) . Eq. (p9)) is suitable for the extrap- 
olation schemes described in Sec. IIIIl since the indices 
G, G' and G" refer again to plane-waves. Inserting a 
masking function that has been introduced for the union 
Eg cutoff into Eq. (pQ]) gives 

i G,G',G" (42) 

xPg(G'-G",G-G";M), 

where 

xfG^'(M') = t^:(2v^ - igfr). (43) 

Note that only three out of four orbital indices are trans- 
formed, and that the transformed x is not symmetric. 
Due to momentum conservation in the transformed basis, 
the (truncated) correlation energies obtained are, how- 
ever, invariant with respect to the transformation of i. 
Note that Ecori-,cs(M' , M) converges towards Eq. (PS)) for 
a sufficiently large M. Replacing Pg(G'-G", G-G"; M) 
with Pfe(G',G";M) = 6 (G') O (G") in Eq. (|42]) yields 
effective basis set energies analogous to the iJ^-cutoff de- 
scribed in Sec. IIIBI 



B. Computational details 

We employ the Vienna ab-initio simulation package 
(VASP) in the framework of the projector-augmented 
wave (PAW) method to carry out MP2 calculations of the 
LiH solid and molecule.— In the PAW method the one- 
electron orbitals ip are derived from the pseudo-orbitals 
ip by means of a linear transformation^^ 

\iJ) = H)+J2m^\^^)){p^\i^). (44) 

i 

The pseudo-orbitals ip are the variational quantities of 
the PAW method, and are expanded in reciprocal space 
using plane waves. The index i is a shorthand for the 
atomic site R^, the angular momentum quantum num- 
bers li and mi, and an additional index denoting the 
linearization energy^. The all-electron partial waves 4>i 
are the solution to the radial Schrodinger equation for the 
non-spin-polarized reference atom at specific energies Ci 
and specific angular momentum li. The pseudo-partial 
waves, </),;, are equivalent to the all-electron partial waves 
outside a core radius Tc and match continuously onto (pi 
inside the core radius. The partial waves (pi and (pi are 
represented on radial logarithmic grids. The projector 
functions pi are constructed in such a way that they are 
dual to the pseudo partial waves, i.e., 

(45) 

For a more detailed outline of the PAW method and a 
thorough discussion of the evaluation of electron repul- 
sion integrals in VASP we refer the reader to Ref J^. 



The employed plane wave basis set for the one-electron 
orbitals and the transformation matrix U is defined by 
all PWs e**^"" with wavevectors G satisfying the equation 



For the calculations of LiH we use E'cut = 400 eV. The 
evaluation of electron repulsion integrals in the PAW 
method requires an auxiliary plane wave basis set. We 
choose our auxiliary plane wave basis set to be identical 
to the basis set defined by E'cut- 

In the present work, we employ 200 and 50 natural 
orbitals to calculate the correlation energies of the solid 
and molecule, respectively. Convergence in the natural 
orbitals basis is two times faster than using Hartree-Fock 
orbitals. Natural orbitals are calculated by diagonalizing 
the one-electron reduced density matrix. A detailed ex- 
planation of this procedure can be found in Ref J^. 

For the LiH solid calculations, we employ a supercell 
containing 8 Li and 8 H atoms. The supercell has a 
volume of 136.24 A"^. The LiH molecule is simulated 
using a box with a volume of 91.12 A'^ and a bond length 
of 1.595 A. The Li Is electrons are frozen and do not 
contribute to the correlation energies. 

C. Results: LiH molecule and solid 

In the following we will apply three different cutoff ex- 
trapolation schemes to the LiH solid and molecule using 
MP2: (i) the local Eg cutoff extrapolation scheme that is 
equivalent to the one previously outlined in Ref. [l^ (ii) 
the union Eg cutoff, and (iii) the iJfc-cutoff. Figures[H]and 
|9] show the convergence of the MP2 correlation energy of 
the LiH solid and molecule, respectively. Both single- 
point extrapolations show a much-improved behaviour 
over the previous scheme that is analogous to a local Eg- 
cutoff, where arcing causes pathological behaviour and 
poor CBS estimates at low M. In both solid and molec- 
ular LiH, the (union, SPE) £'g-cutoff seems to converge 
quicker. 

VI. CONCLUDING REMARKS 

In this paper, we have investigated the convergence 
of correlation energies using plane-wave wave-function 
expansions. Starting by treating the finite simulation- 
cell electron gas with the simplest correlated quantum 
chemical method, second-order M0ller-Plesset perturba- 
tion theory, we derive a functional form of the finite basis 
set correlation energy of 1/M, where M is the number 
of plane waves enclosed by a spherical cutoff in fc-space. 
Although perturbation theory diverges in metallic sys- 
tems for any strength of Coulomb interaction, the quali- 
tative behaviour of the wavefunction around the correla- 
tion hole is in common with other higher-level methods. 
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FIG. 8: The MP2 correlation energy of the LiH 2x2x2 
supcrccU retrieved as a function of the basis set size for 

a variety of extrapolation schemes. The SPE curve 
shows the effective basis set energies produced from a 
single calculation with M' = 2045. 



the complete basis set limit correlation energy in CCD 
and RPA+SOSEX. 

Finally, this is applied to real materials (molecular and 
solid LiH). We find that the energies computed by single- 
point extrapolation converge better and more reliably 
than previous extrapolation techniques'^. It is our hope 
that this can be applied in future plane-wave wavefunc- 
tion based calculations. 
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Appendix A: CCD and RPA+SOSEX 

In the following, we will briefly outline coupled- 
cluster doubles (CCD) theory and the random-phase 
approximation plus second-order screened exchange 
(RPA+SOSEX). 

CCD is a widely used quantum chemical method to 
study the electronic ground state energy of atoms and 
molecules and relies on an exponantial Ansatz for the 
many-electron wavefunction that reads^L^ 



CCD 



(Al) 



FIG. 9: The MP2 correlation energy of the LiH 
molecule retrieved as a function of the basis set size for 

a variety of extrapolation schemes. The SPE curve 
shows the effective basis set energies produced from a 
single calculation with M' = 1647. 



We verify that this behaviour extends to coupled- 
cluster doubles (CCD) and the random-phase approxima- 
tion plus second-order screen exchange (RPA+SOSEX), 
in common with exact results from full configuration in- 
teraction quantum Monte Carlo (FCIQMC)^>2^. 

By viewing the distribution of the wavefunction in con- 
figuration space over double-excitations, and relating this 
to orbital momenta in fc-space, we propose several new 
basis set truncations based on the momentum transfer 
vector. We discuss these in terms of their comparative 
speed and smoothness of convergence recommending one 
scheme, which we call the union Eg-cntoS, that gives 
overall the most desirable properties. This is then gener- 
alised to other single-reference quantum chemical tech- 
niques, allowing for the development of a single-point 
extrapolation technique which uses information from a 
single large-basis-set calculation to provide estimates for 



where T2 refers to the double excitation operator.— 

occ. 

T2|*o) =^ 



OCC. unocc. 

i<j a<b 



occ. unocc. 
i<j a<b 

(A2) 

We choose to be the Hartree-Fock reference deter- 
minant. The solution to the CCD wavefunction is ob- 
tained by projecting vj/^:^^ onto a set of doubly excited 
determinants. This set of equations is termed amplitude 
equations. The CCD amplitude equations read^ 



0=< + (C + 

+E[(2<c'- 



,HF 



,HF 



aUbc 
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bl ^ac 
IV -Lab 



bUac 
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and a, b, c refer to occupied and unoccupied orbitals, 
respectively. The amplitude equations can also be writ- 
ten in a more compact fashion by defining intermediate 
quantitiesi^. Solving Eq. (|A3[) yields the wavefunction 
coefficients in configuration space tf^ and allows for the 
correlation energy to be calculated according to Eq. ([35|l . 
Due to the computational cost involved, CCD has so far 
only rarely been applied to solid state systems. 

Freeman, and Bishop and Liihrmann studied the 
uniform electron gas using an approximation to CCD 
theory ^^1'^^ This approximation has recently attracted 
renewed interest and is termed RPA+SOSEX,^"i^ 
RPA+SOSEX differs from CCD in two points: (i) the HP 
reference is replaced by the KS reference, which greatly 
reduces the one-electron gap and, (ii) the double ampli- 



tude equations are approximated by so-called ring dia- 
grams only 

+ T.''?ct^+Y.*>^ + Y.^>c:4': (A4) 

Ic Ic ll'cc' 

Once obtained, the RPA-hSOSEX tfj'-amplitudes can be 
employed to calculate the RPA-I-SOSEX correlation en- 
ergy using Eq. (|35p . A rigorous justification for this ap- 
proximation is not straightforward and would be beyond 
the scope of this work. However, Ref. IH outlines the 
connection between the above amplitude and Casida's 
equation. 
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